Continuous RSB mean-field solution of the Potts glass 
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We investigate the p-state mean-field model of the Potts glass (2 < p < 4) below the continuous 
phase transition to a glassy phase. We find that apart from a solution with a first hierarchical 
level of replica-symmetry breaking (IRSB), locally stable close to the transition point, there is a 
continuous full replica-symmetry breaking (FRSB) solution. The latter is marginally stable and has 
a higher free energy than the former. We argue that the true equilibrium is reached only by FRSB, 
being globally thermodynamically homogeneous, whereas IRSB is only locally homogeneous. 

PACS numbers: 64.60Cn,75.50Lk 



The Potts model is a generalization of the Ising model 
to more than two spin components. It has been a sub- 
ject to intense research in recent decades. The origi- 
nal formulation of Potts [T] with a Hamiltonian Hp = 
~ Ji]^n,,nj where = 0, . . . ,p - 1 is an admissi- 

ble value of the p-component model on the lattice site 
Ri, is unsuitable for practical calculations. The Potts 
Hamiltonian can be represented via interacting spins [2] 

H.p = --'^Jij^i-^j , (1) 

where = {s\, . . . sf~^} are Potts vector variables tak- 
ing values from a set of state vectors {^aW^i- Functions 
on vectors are in equilibrium fully defined through 
their scalar product • = 6^46^ — pSab — 1, 
a € {1, ...,p—l}. We use the Einstein summation conven- 
tion for repeating Greek indices of the vector components 
indicating a scalar product of Potts vectors. 

Interest in the Potts model was recently revived by 
discovering its interesting properties when frustrated 
randomness is introduced. If we assume that the ex- 
change Jij is a random variable we speak about a Potts 
glass. Its mean-field (long-range) version with a Gaus- 
sian distribution P{Jij) = (27r J^/iV)"^/^ exp{-( J,y - 
Jq/NY /2,P /N} is a generalization of the Sherrington- 
Kirkpatrick (SK) model of a spin glass [3]. Analogously 
to the pure model, the Potts glass shows a richer crit- 
ical behavior than the Ising spin glass (p = 2) [U [5]. 
The peculiarity of the mean-field solution of the p-state 
Potts glass for p > 2 is that the replica symmetric (RS) 
solution below the transition point to the glassy phase 
does not break into the Parisi solution with a continu- 
ous replica symmetry-breaking but rather to a first RSB 
hierarchical state [6, 7J. The existence of a locally sta- 
ble IRSB solution generated from an instability of the 
replica-symmetric one may enlighten the construction of 
the the equilibrium glassy phase. In the effort to solve 
the SK model, Parisi proposed a hierarchical construction 
how to break the replica symmetry when using the replica 
trick to average over the random spin couplings 8J. Al- 
though it is always necessary to construct the equilibrium 



glassy state via discrete (integer) numbers of replicas, so- 
lutions with finite many replica hierarchies were for long 
assumed as auxiliary, intermediate states. Only the con- 
tinuous Parisi solution was considered to represent the 
true equilibrium. This was the case in the SK model, 
where all iCRSB solutions proved unstable. Local stabil- 
ity of the IRSB solution of the Potts glass renewed in- 
terest in solutions with finite many hierarchies of replica 
generations. 

The SK model proved very inspiring in understanding 
the low-temperature glassy phase. Not only the replica- 
symmetry breaking of Parisi was derived within it but 
also a mathematical poof was presented that the con- 
struction of Parisi leads in this model to an exact equi- 
librium state [HI Hn]- However, the Ising spin glass has 
not allowed for understaning of the actual meaning of the 
states with finite many hierarchies of replica generations. 
The rigorous construction of Guerra and Talagrand does 
not specify how many replica hierarchies are necessary to 
reach the maximum of free energy. On the other hand, 
the p-state Potts model of a glass can be used to clar- 
ify this point as well as the way how and when, if ever, 
iCRSB solutions break into the Parisi FRSB one. It is by 
now clear that the IRSB solution of the Potts glass does 
not stay stable at all temperatures [6j[12]. It breaks down 
at a lower temperature Tq to another state as discussed 
in the framework of the p-spin glass [13] . It is generally 
assumed that the equilibrium state below this tempera- 
ture is the Parisi FRSB solution, though not rigorously 
proven. Namely, it remains unclear when and how the 
FRSB solution emerges in the p-state Potts glass or the 
p-spin glass. 

It is the aim of this Letter to trace the genesis of 
the continuous replica-symmetry breaking solution of the 
mean-field p-state Potts glass. We show that it emerges 
at the transition point to the glassy phase and coexists 
with the IRSB solution. Although the former is only 
marginally stable and the latter locally stable, the FRSB 
solution has a higher free energy and unlike IRSB it is 
globally thermodynamically homogeneous. We demon- 
strate by explicit calculations near the transition point 
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to the glassy phase that IRSB is globally unstable and 
FRSB offers the true equilibrium with the extremal free 
energy when the replica-symmetric solution becomes un- 
stable. 

Thermodynamic homogeneity was introduced in spin 
glass models by one of us [H] as a condition to be fulfilled 
by free energy in the thermodynamic limit. Instability 
of the replica-symmetric solution and replica-symmetry 
breaking in the replica trick were shown to be equivalent 
to a breakdown of thermodynamic homogeneity. Thereby 
a physical interpretation was given to replica-symmetry 
breaking. When the statistical system is thermodynami- 
cally inhomogeneous the thermodynamic limit of an iso- 
lated system (micro-canonical ensemble) differs from the 
thermodynamic limit of a system embedded in a ther- 
mal bath (canonical ensemble). The thermal bath can 
be generated by scaling or replicating the original phase 
space. We probe a dependence of the original system on 
the bath variables via an infinitesimal interaction with 
the replicated ones. Thermodynamically homogeneous 
systems must not be macroscopically affected by such 
perturbations in the thermodynamic limit. If the original 
system is macroscopically influenced by an infinitesimal 
interaction with the bath, new thermodynamic variables 
(chemical potentials) must be introduced to equilibrate 
the system and the bath [15] . Local thermodynamic ho- 
mogeneity guarantees only independence with respect to 
small (infinitesimal) scalings of the phase space, that is, 
situations with the bath in a thermal equilibrium. We 
show that the locally stable IRSB solution appears glob- 
ally unstable and depends on the variables of higher RSB 
solutions if the higher spin hierarchies (bath) are pulled 
out of equilibrium. Only global thermodynamic homo- 
geneity, satisfied by FRSB, warrants that the equilibrium 
state does not depend on the behavior of the bath (scaled 
or replicated variables) in which the thermodynamic sys- 
tem is embedded, or with which it can interact. 

The starting point for the construction of the equilib- 
rium state for spin-glass systems is a free energy with K 
spin hierarchies. We will deal with the mean-field (long- 
range) p-state Potts model with Hamiltonian Hp from 
Eq. ([ij and neglect all spatial fluctuations. We assume 
a Gaussian randomness in the exchange parameters Jij 
with Jo = for p < 4 so that the transition to the glassy 
phase is higher than that to a ferromagnetic state [4]. 
A hierarchical free energy with K hierarchies depends 
on 2K -\- 1 variational parameters. It is a parameter q 
weighing the fluctuations of the random exchange and K 
pairs Ax/ , mi governing the interaction between the hi- 
erarchies of replicated spin variables [TJ . Mean-field free 
energy with K hierarchies for the p-state Potts glass has 



the following form 
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J Ppylog4(y) (2a) 
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The initial zero-level partial sum reads 
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We denoted a vector Gaussian 



exp<(/?J ( V^y"-^^v^A° 

(2c) 

measure 'Dpy = 

Ur=\dyo.e^p[-iy»f /2} /V2^ and {A}, 

{A;, A;+i, . . . Xk}- Notice that each integration variable 
y, Xi is a p — 1 dimensional vector. It is worth noting that 
partition sum zq depends on a particular representation 
of the basis vectors e/, but the physical quantities, 
averaged over fluctuating vector fields y and A;, are 
invariant with respect to gauge transformations in the 
representation space of the Potts spin-like variables. 

The rigorous construction of Guerra and Talagrand 
tells us that the equilibrium state is reached by supre- 
mum of functionals fx from Eq. ^ among all possi- 
ble choices of parameters K] q, Axi,mi. It, however, 
does not say whether the number of hierarchies is finite 
or must be infinite. The Parisi continuous symmetry- 
breaking solution within this construction is obtained 
as an asymptotic limit K ^ oo with an assumption 
of a homogeneous distribution of overlap susceptibilities 
Axi 1 /K [T5] . Each hierarchical free energy satis- 
fying K stability conditions represents a local maximum 
if < TO/ < 1 [14]. The question is whether a locally 
stable solution with a finite number of hierarchies is also 
the global maximum. 

One of us recently derived a representation of the 
Parisi free energy with FRSB in closed form being decou- 
pled from the iterative construction via discrete replica- 
symmetry breaking hierarchies [17]. This representation 
allows us to look for a solution with FRSB without inves- 
tigating local stability of the discrete iCRSB solutions. It 
is a straightforward task to generalize the representation 
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from Ref. TTto the p-state Potts glass. We obtain 

- [q, X, to(A)] ^\ogp+^{p-l){l-q-Xf 

-^{p-l)X J' dXm{X)[q+X{l-X)] + {g{l, h + y^))^ 

(3a) 

where 

g{iy, h) = Ta exp I y ^ dX [dj,^ dj,o, 

+miX)g'^iX,h + h)dj,^] |go(h + h)|^^„ (3b) 

and 5o(h) = lnX]f=i exp{/3/i"e"} is the local interacting 
part of the Potts free energy. We introduced an evolution 
operator represented via a "time-ordering" operator 
ordering products of A-dependent non-commuting oper- 
ators from left to right in a A-decreasing succession. We 
further denoted 5^(A, h) = dh^g{X, h). 

Notice that the FRSB state is represented only by a 
single function, m(A), A G [0, 1] in our case. It generalizes 
the sequence mi from the discrete scheme. It is related 
to the Parisi order-parameter function via m{X) = x{q + 
X{1 — X)). The overlap susceptibilities Axi transform to 
XdX. The Edwards- Anderson parameter here is qEA = 
q+X. That is, our function m(A) is essentially the inverse 
Parisi order-parameter function (with the inverse slope). 



We now have two free energies, Jk from Eq. ([2| and 
fc from Eq. ([3]), for which we can try to find stationary 
solutions in the space of their variational parameters. It 
is clear that we are able to resolve these free energies and 
the stationary order parameters only near the continu- 
ous transition to the glassy phase. To find the transition 
point we start with the RS solution with the only pa- 
rameter g j4j . If we denote the transition temperature Tc 
and T = (Tc — T)/Tc we obtain in the leading order of t 
the asymptotics q^s = 4r/(6 — p). We chose the energy 
scale J = 1. The stability condition for the RS solution 
is Ao = 2t(2-p)/(6-p) [7]. The RS solution is unstable 
(Ao < 0) in the leading order of r for p > 2. The RS 
solution of the Ising spin glass {p = 2) gets unstable only 
in t2 [TB]. 

We asymptotically expand free energies of the IRSB 
and FRSB solutions in powers of the small parameter r. 
We will calculate only the differences of the free energies 
to the paramagnetic one, being fSfpara — — /3^/4(p— 1) — 
log(p). It is a rather tedious task to do so and the details 
of this calculation will be presented elsewhere. We only 
note here that we used the program MATHEMATICA 
to generate the expansion coefficients of the solutions for 
the order parameters. We had to expand them up to the 
second order for IRSB and to the third order for FRSB 
so that to obtain the expansion of the free energies to the 
fifth order where the two solutions differ. 

For IRSB we derived the following asymptotic free en- 
ergy 
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The variational parameters for IRSB in the leading 
asymptotic order are q = and Ax == 2t/(4 — p), 
mi = (p-2)/2-h(36-12p-|)2)r/4(4-p). The transition 
from the RS to the IRSB solution becomes discontinu- 
ous for p > 4. Unlike the RS solution, the IRSB solution 
becomes locally stable for p > p* as can be deduced from 
the stability criterion explicitly evaluated near the tran- 
sition point 



Ai = 6(43^(^-24^+12) 



(5) 



The stability function Ai becomes positive for p > p* ^ 
2.82. 

The expansion of the FRSB solution below the tran- 
sition point to the glassy phase is more intricate due 
to a complicated structure of representation in Eq. (|3]). 
We have to solve asymptotically equations for the order- 
parameter function m(A) for each A G [0, 1]. We find that 
q = identically and X = 2t/(4— p) recovers in the lead- 



ing order Ax from IRSB. Finally the order-parameter 
function in the leading-order reads 

m(A)-^ + i[(3-2p)p + 6]A 



1 



i7p - 2A)p + 12] XX . (6) 



Expansion of the free energy with the appropriate asymp- 
totic solutions for the order parameters then is 
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(7) 



Comparing the asymptotic free energies for the IRSB 
and FRSB solutions we obtain a difference 



P{fc ~ IiRSb) — 



(p-l)(p(7p-24) + 12) 
720(4 - pf 
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(8) 
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FIG. 1. DifFerences in the leading asymptotic order below the 
transition point to the glassy phase between the free energies 
of the FRSB and IRSB (solid hne) and FRSB and RS (dashed 
line) solution. Notice that the FRSB solution splits from the 
RS one at the order r'^ (right scale) and form IRSB not till 
the fifth order (left scale). 



telling us that the free energy from FRSB has for the 
Potts glass a higher free energy than that from IRSB for 
2 < p < 4. Potts glass for p > 4 demands a more detailed 
analysis, since the denominator on the right-hand side 
of Eq. ^ goes through zero. We also calculated the 
difference of the FRSB solution to the RS one leading to 



/3(/c - fBs) 



2^3 



3(4-p)(6-p)2 



(9) 



Notice that /c = fiRSB to this order. We see again that 
our asymptotic expansion is meaningful only for p < 4 
where the transition to the RSB phase is continuous and 
T serves as a small parameter. The 4-state Potts glass 
is a singular (degenerate) case, since mi = 1 at of 
the IRSB solution. Both differences in free energies are 
plotted as a function of p in Fig. [T] It is interesting to 
notice that /i = fc for p — p* above which the IRSB 
solution becomes locally stable. The free energies equal, 
however, at this point only in the leading order and this 
degeneracy is lifted in a higher order of r. 

Last but not least, the FRSB solution of the Potts glass 
suggests that the physical interpretation of the Parisi 
order-parameter function q{x) should be reconsidered. 
Parisi [T5] proposed that the derivative dx{q)/dq be in- 
terpreted as the average of the probability distribution 
P{q) of different pure states with overlap q. This inter- 
pretation seems not to work for the Potts glass, since 
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where q = X{1 — A). Function P{q) gets negative for 
p > p* K, 2.82, the stability region of the IRSB solution. 

To conclude, our analysis of the asymptotic solution of 
the p-state Potts glass for 2 < p < 4 below the continuous 
transition to the glassy phase revealed, against general 
expectations, that the IRSB solution is not the true equi- 
librium, although locally stable for 4 > p > p* « 2.82. 
We found a different state with FRSB that continuously 
develops from the critical point T^. The FRSB solution 
has a higher free energy than that of the IRSB one and 
represents the global maximum as demanded by the rig- 
orous construction of Guerra and Talagrand. Below the 
critical point T < for 2 < p < 4, the IRSB solu- 
tion is globally unstable and decays on a long but finite 
time scale to the FRSB solution. We could give the de- 
rived result a physically appealing interpretation in that 
the IRSB solution is only locally thermodynamically ho- 
mogeneous while the necessary condition of global ther- 
modynamic homogeneity is fulfilled only by the solution 
with the continuous FRSB. 

Research on this problem was carried out within 
project AV0Z10100520 of the Academy of Sciences of the 
Czech Republic. 
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